Reduction of surface coverage of finite systems due to geometrical steps 
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The coverage of vicinal, stepped surfaces with molecules is simulated with the help of a two- 
dimensional Ising model including local distortions and an Ehrlich-Schwoebel barrier term at the 
steps. An effective two-spin model is capable to describe the main properties of this distorted Ising 
model. It is employed to analyze the behavior of the system close to the critical points. Within 
a well-defined regime of bonding strengths and Ehrlich-Schwoebel barriers we find a reduction of 
coverage (magnetization) at low temperatures due to the presence of the surface step. This results 
in a second, low-temperature transition besides the standard Ising order-disorder transition. The 
additional transition is characterized by a divergence of the susceptibility as a finite-size effect. Due 
to the surface step the mean-field specific heat diverges with a power law. 

PACS numbers: 64.60.Fr, 64.70.Nd, 75.70.Ak, 68.35.Rh 



I. INTRODUCTION 

The characterization of phase transitions becomes es- 
pecially demanding in situations where the order param- 
eter is not directly accessible by experiment. One ex- 
ample is the search for a nuclear liquid-gas phase transi- 
tion. A considerable discussion can be found in the lit- 
erature about the possibility to observe a negative heat 
capacity as one signal of such a possible liquid-gas phase 
transition! Such negative heat capacities appear in finite 
systems which are adequately described within the micro- 
canonical ensemble. We report here an observation that 
a transition with a divergent heat capacity can occur as a 
consequence of geometrical distortion rather than due to 
a phase transition even in a canonical treatment. This 
may shed some light on the nature of observed signals. 

One meets a similar situation when describing the cov- 
erage of surfaces with molecules. There it is interest- 
ing to distinguish signals caused by phase transitions 
between different adsorbate arrangements from signals 
due to structural transitions at local deviations from the 
ideal surface geometry. Different surface defects have 
been studied within Ising models^ by density renormal- 
ization methods as well as Monte Carlo techniques. Non- 
universal features were observed and the critical expo- 
nent of the magnetization was found to be near 1 /2 for 
infinite systems. A review on the vast literature about 
phase transitions in inhomogeneous systems can be found 
ml 

We investigate here a finite-size two-dimensional Ising 
model suitable to simulate the coverage of surfaces by 
molecules. While the explicit simulation with realistic 
parameters was described in jj| we concentrate here on 
principal results how the surface modification is influenc- 
ing the transitions and the critical exponents. We sug- 
gest that the occurrence of divergent (or negative) heat 
capacities is not a unique signal of a phase transition but 
can occur due to the geometrical distortion of the system 
accompanied by anomalous exponents, which even fulfills 



the scaling hypothesis. 

The two-dimensional Ising model belongs to the most 
studied models. For an overview see |6(]. The exact 
solution 7 shows a phase transition with a critical behav- 
ior: 

spontaneous magnetization M ~ |T — T c \@ 
magnetic field dependence H ~ \Mt=t c \ S 
susceptibility \ ~ \T - T c \-~> 

specific heat c H ~ \T - T C |~ Q . 

Two exponents are exactly known, i.e. [3 = 1/8 and 
7 = if 0- From asymptotic expansions and strong 
numerical evidence one has furthermore a = and 
5 = 15 01 where the specific heat diverges logarithmi- 
cally Weiss' mean-field approximation instead leads to 
a = 0, /3 = |, 7 = 1, and <5 = 3 [Till. Both sets of critical 
exponents fulfill the inequalities 12] a + 2/3 + 7 > 2 and 
[l3j a + (3(1 + 6) > 2 known as scaling hypothesis. These 
scalings are determining the corresponding universality 
classes with specific scaling functions for the magnetic 
field dependence of the magnetizationiS*i^. The univer- 
sality of this phase transition in two dimensions has been 
experimentally confirmed 15 . Recently, the universality 
has been investigated with respect to finite size scaling 16 
and oscillating fieldsit 

Modifications of the scaling relation due to surface de- 
fects have been studied extensively, see citations in [l8| . 
Let us only mention some of the results found in print. 
The divergences of the specific heat for free and ferromag- 
netic boundaries in different Ising models have been stud- 
ied for 40 years 2 * The effect of a surface in an Ising model 
induces spatial correlations which could be treated with 
the help of a Ginzburg-Landau equation^. The two-spin 
correlations induced by a line defect in an Ising square 
lattice are considered with the help of two-particle cor- 
relation functions^* 2 ! Many-point correlation functions 
along a modified bond have been calculated as well 22 . 
The critical exponents for the magnetization of a line 
defect are known analytically 2 ^^ 5 -. 

We will present here a quadratic Ising model with a 
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line defect and an additional change of the magnetic 
field along the line known as Ehrlich-Schwoebel barrier. 
This can mimic the surface coverage with molecules in 
the presence of an additional step. First we explain the 
model and present the numerical results. We will find 
that the Ehrlich-Schwoebel barrier induces an additional 
transition. Then in the third chapter we show that the 
commonly used standard mean-field model fails to ex- 
plain the observations quantitatively. An effective model 
is suggested which accounts for the basic results. This 
effective model is then discussed in chapter IV with re- 
spect to the critical exponents and is compared to the 
mean-field exponents of the standard Ising phase transi- 
tion. 



II. ISING MODEL WITH SURFACE STEP 

In order to simulate the coverage of surfaces by 
molecules we imagine this surface as an N x N square lat- 
tice with a straight step across the middle of the lattice. 
The spin-up states describe a molecule sticking to the 
surface while the spin-down states describe the absence 
of a bound molecule. The interaction with the j • = 4 
neighboring surface molecules is described by the cou- 
pling constant J. Across the surface step we choose a 
different coupling constant J' = J/v2. The interaction 
of the surface molecules with the substrate background 
is modeled in analogy to the spin coupling with an exter- 
nal field. Therefore we shall use the external magnetic 
field as a synonym for the coupling of molecules with the 
background. At the sites adjacent to the step the mag- 
netic field is augmented by an additional term, H s > 0. 
H s models the Ehrlich-Schwoebel barrier, which impedes 
the diffusion of adsorbates across surface steps. On top 
of the step edge H s is added to H, hence it locally re- 
duces the attractive adsorbate-substrate interaction and 
mimics the lower density of coordination sites on top of 
the step edge. From below, H s is subtracted, thus it 
enhances the adsorbate-substrate interaction and models 
the higher number of coordination sites along the step. 
Motivated by the results of Ehrlich and Schwoebel on the 
stability of step arrays we chose the attractive and repul- 
sive parts of the barrier equally high. Thus, the Ising 
Hamiltonian for the stepped square lattice reads: 



derived 23 to be 



/?h=o = 7-^-arccos 2 



— tanh 



2(J- J') 



(2) 



This shows that the critical exponent becomes dependent 
on the coupling strength and is therefore non-universal. 

We solve the two-dimensional Ising model with the 
standard Metropolis scheme and the magnetization is 
now employed as measure for the surface coverage with 
molecules, plotted in Fig. ^ One sees that with increas- 
ing external field (or coupling of molecules to the back- 
ground) a smearing of the standard Ising phase transition 
is obtained, which results in high temperature tails. This 
effect is well studied and experimentally confirmed^. 
The ferromagnetic transition occurs only for vanishing 
magnetic fields. Fig. ^ also compares the solution of the 
two-dimensional Ising model with and without a surface 
step: For low values of the external magnetic field the 
step leads to a characteristic reduction of the magnetiza- 
tion at low temperatures. 




FIG. 1: Magnetization versus temperature scaled by the half- 
width temperature where the magnetization is 1/2 for differ- 
ent external magnetic fields. The Ising lattice (50 x 50) with 
a step (open symbols and dashed lines) is compared with the 
one without a step (closed symbols and solid lines). The pa- 
rameters are J = 1, J' = l/v2 and fioH s — 5. Error bars are 
of or less than the size of the symbols. 



H = JsiSj + J'si'Sj' + Js 



i',3 



III. CONSTRUCTION OF AN EFFECTIVE 
MEAN-FIELD MODEL 

A. Local mean-field model 



where i sums over all N(N — 2) terrace sites, i' over the 
step sites, j over all neighbors with coupling J and j 
over all neighbors with coupling J'. Without fields H 
and H s this Hamiltonian is an Ising square lattice with a 
ladder defect^ and the exact critical exponent has been 



This behavior can be understood by an effective two- 
spin model. We describe briefly in the following that the 
standard mean-field model, as e.g. used in the appendix 
of yi|! fails. For a lattice size of N x N spins and pe- 
riodic boundary conditions, the system is homogeneous 
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in the direction parallel to the step, thus we can restrict 
our considerations to the direction perpendicular to the 
step. We distinguish N — 2 normal Ising spins and 1 spin 
at each side of the step, for which the modified coupling 
constant J' and the Ehrlich-Schwoebel barrier term H a 
are taken into account 26 . All three kinds of spin experi- 
ence an effective mean-field. We will denote the normal 
mean spin with s and the mean spins at the step with 
according to the sign of the Ehrlich-Schwoebel barrier 
term, ±H S . N — 4 of the sites occupied by normal spins 
see a mean-field H consisting of the external field H and 
the interaction with 4 neighbors, 4J. The remaining two 
of the N — 2 normal spins interact only with 3 normal 
spins and with one spin at the step. Therefore we have 



JV-4 

LIqH = (j, Q H + — - 2 4Js 



2 qT -L q~ 

- (3 J S + J^^).( 3 ) 



N-2 



The spins along the step have two interactions with the 
same kind of spins, 2Js ± , one neighbor with normal 
coupling, Js, a contribution from the coupling across 
the step, J' , and an interaction with the substrate of 
H ± He. This results in 



l^o^ = fi H ± fi H s + 2 Js ± + J's T + Js. 



(4) 



The partition function is then trivially written as 

z = [2cosh(/3^ 0J ff)] (Ar ~ 2)Ar 

x [2 cosh (/3fi H + )] N [2 cosh (f3^i H-)} N (5) 

with the inverse temperature (3 = 1/ksT. The mean 
spins are calculated by expressions of the statistical av- 
erages, s — ta,nh(Pfi H) and = tanh (/3fi H ± ). © 
and (@J represent the self-consistent mean-field equations 
for s and s . This mean- field result is exactly equivalent 
to the Bragg- Williams method by minimizing the Gibbs 
functional and assuming that the many-spin correlation 
function factorizes into single-spin ones. Such mean-field 
equations for open surface defects have been investigated 
in [H|l|. 

First it is instructive to solve this equation in the limit 
of zero temperature. Then one gets the values of the 
mean spins ±1 according to the sign of the mean fields 
© and Consequently the total mean spin 



M 



1 



2 

N 



1 



+ tt( s +s ) 



N 



(6) 



1 



f for T -> if 



approaches the reduced value M 
H > and (H+ ^ Oand/J- ^ 0) . Therefore, as seen 
in Fig. ^ the reduction is 1 — 2/50 = 0.96 at low temper- 
atures and low external fields. Discussing the different 
cases and taking into account that the partition function 
assumes the maximum one deduces that this reduction 
happens if and only if 



fj,oH + 3 J + J' < hqH s 



(7) 



as outlined in the appendix. Though this mean-field 
model obviously describes the reduction qualitatively the 



actual numbers do not agree with the simulation result. 
Therefore we can conclude that the standard mean-field 
model is not able to describe the effect quantitatively. 
This is understandable since surface defects induce non- 
local correlations^. These nonlocal correlations result in 
a spatial dependence of the magnetization on the dis- 
tance from the step on the surface. This can be modeled 
by a Ginzburg-Landau equation as derived in jl9( ■ 

B. Effective mean- field model 

A better match with the numerical data is obtained 
for an effective two-spin model taking into account these 
nonlocal correlations in a certain sense. We discriminate 
now only between normal spins s on attractive sites and 
fictitious spins s' at the repulsive sites with H — H s along 
the step. In this approach, each row across the terrace 
contains N — 1 sites with normal spins s and the repul- 
sive site with spin s' . An energy-conserving mapping of 
the intuitive three-spin model described above onto this 
simplified two-spin model is possible by setting ,s + = s 



and s = 



This mapping relies on the following con- 



siderations: The presence of one and only one normal 
spin type s is only guaranteed if the effective field is ho- 
mogeneous on the terrace sites. On the other hand, the 
two-spin model explicitly accounts only for the repulsive 
part of the Ehrlich-Schwoebel barrier along the step, and 
omits the attractive part. Yet, the attractive part of the 
Ehrlich-Schwoebel barrier must not be neglected, as it is 
employed to obtain the data from the numerical simu- 
lations. The only solution consistent with both require- 
ments is to distribute H s evenly over the terrace as an 
overall enhancing field of the strength H S /(N — 1). When 
summing over all N — 1 terrace sites in the Hamiltonian, 
the same total energy of the system results for the two- 
spin and the three-spin models. 

The condition for the second, low-temperature transi- 
tion can then be obtained by calculating the site energy 
E 1 of the position with spin s' as total energy difference 
between the three-spin model and the attractive part of 
the two-spin model. This procedure yields for the site 
energy 



E' = fi H{s 
(* + ) 2 



-2J 



+ (s-y + 



s) + n H s (s + - s - s) 
1 



ss^ + ss 



2s z 



+2JV 



(8) 



Equating s — s + — 1 and 



for the site energy E' 



- s~ = — 1 one obtains 
-HqH + {.iqH s — J — J' for the 
effective spin s' in anti-parallel orientation with respect 
to s. This orientation is favorable, if E' < 0, hence the 
condition for a reduction of the magnetization reads: 



li H + J + J' < hqH s 



(9) 



Indeed, numerical simulations for different parameter 
sets J, H, H s confirm this result. Hence, the effective 
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two-spin model is employed for the further analysis of 
the numerical simulations. 

The mean field of the normal spins is calculated anal- 
ogously to © 



HoH = /i H 



TV- 3 

: L 

N — 1 



-AJs 



6s + 2s' 
N-l 



J = [m,H + j Js(10) 



The fictitious spins differ from the normal ones by a 
constant s' = c + s where — 2 < c < 2. For a given 
value of T, each spin s has therefore a mean field jJs = 
AJs + 2Jc/(N — 1) sa AJs with a maximal relative error 
of 1/{N — 1) which we can neglect in the following. 

The effective spins along the step are described by a 
mean field consisting of the linear combination of the 
couplings J, J' with the spins s,s'. Taking into account 
© and that a possible second transition can only occur 
at a temperature T' r ~ J' as well as that for J = J' 
and H s = the normal Ising model should reappear, we 
obtain uniquely the mean field of the effective spins as 



J') 



2JV 



fill 



fi H' = hqH - fioH s + (J 

the derivation of which is outlined in the appendix. The 
partition function can again be trivially written and the 
mean spins are 

s' = taxLtL\j3(fio(H-H s ) + (J + J')s + 2J's')} 
s = tanh [/3(/i i? + jJs)], (12) 

representing the self-consistent mean-field equations for 
s and s'. We find from l|12|l that the total mean spin 



M = 



1-1 

N 



1 , 

S+ N S 



(13) 



approaches the reduced value M = 1 — 2/iVin the limit 
of zero temperature if and only if the condition (JHJ is 
fulfilled as shown in the appendix. 

The solution of Eq. i|12|) versus temperature is plot- 
ted in Fig. [21 for different values of H' . We see the 
characteristic reduction of the effective magnetization 
to 1 — 2/10 = 0.8. This reduction occurs as long as 
0.1 + 0.071 = 0.171 < fi Q H s according to ©. Fig. 0c) 
displays that slightly below the critical value {^qH 8 = 
0.17) we do not have a reduction at T = but a sharp 
drop of the magnetization around T/T c = 2J'/j J = 0.35. 
This is related to a maximum in the specific heat at 
a second critical temperature T' c 0.35T C besides the 
usual Ising transition temperature T c as shown in Fig. [31 
The same second transition appears in Fig. ^ where 
HoH s = 5 and consequently the reduction occurs as long 
as hqH < 3.29. 



IV. MEAN-FIELD CRITICAL EXPONENTS OF 
NEW TRANSITION 

We can understand the second transition by expanding 
(|12fl for low fields H and calculating the susceptibility 

X=— |i, =0 . (14) 




FIG. 2: The mean spins s (dashed line) and s (dotted 
line) as solution of 1121 versus temperature together with 
the total mean magnetization, M s t cp , (solid line) for fioH s = 
0.10,0.16,0.17,0.18 (a-d). The magnetization without step 
agrees with s. The lattice size is N = 10, the couplings 
J = 0.1 and J' = J/V2, and the external field H = 0. 



From (|12[) we obtain 



• s "> ^ + Tc §^ 



(15) 



which is easily solved and employed to calculate l|14fl . We 
discuss this susceptibility explicitly near the two transi- 
tions. At the usual phase transition T c = jJ/ks where 
s = we have 



X| s =o = 



lim y| 




1 i_ s '2T+iT c '+(i-l)T c \ ^/ kl 



N 



T-{l-s l2 )Tt 



T — T r . 



T-T c 
(16) 



and the typical critical exponent 7 = 1 occurs for finite 
and infinite lattices. The step of the surface does not 
change the critical scaling of the susceptibility. 

Near the second transition at T' c — 2J'/k B where s' = 
we obtain 



Mo 
k B 



l-N 



1-s 2 



N T-(l-s 2 )T c 
1 T+(l- S 2 )((l/j-l)T c 



2 c) 



N 



(T-(1- S 2 )T C )(T-T L 



lim xU'=o 



Mo(l - s 2 )/k B 
T-(l-s 2 )T c 



(17) 



dH 



and for a finite lattice (N < oo) we see that x ~ 
1/(T— T' c ). Consequently, at the second critical temper- 
ature, T' the susceptibility diverges and a sharp drop of 



5 




5 10 

T [k B /J] 



15 



FIG. 3: The specific heat for parameters of Fig. for the 
numerical data of the Ising model with a step. 



magnetization occurs with the critical exponent 7' = 1. 
This second critical temperature does not appear for in- 
finite lattices since the term with 1/(T — T' c ) vanishes in 
the limit N — > 00. We hence conclude that the second 
transition occurs due to the finite spacing between two 
adjacent surface steps. 

Even the quantitative behavior of the mean-field model 
agrees remarkably well with the numerical solution if we 
scale to the corresponding half-width temperatures as 
can be seen in Fig. Especially the low-temperature 
behavior and the drop at the second, low-temperature 
transition at T' c are well described. Since the mean-field 
approximation does not yield the correct critical expo- 
nents of the standard order-disorder transition of the pla- 
nar two-dimensional Ising model it is in accordance with 
previous findings that deviations occur for temperatures 
higher than Tyi<x- 

In order to substantiate the picture of a second, low- 
temperature transition we investigate the remaining crit- 
ical exponents. We find the magnetic field dependence of 
the magnetization by rewriting (|12|l 



k B T 



Tc 
T ' 



1 



■In 



1 



and near T c 
k B T c 



N 



N - 1 



-M 



N - 1 



(18) 



+ o(s 5 ) (19) 



where lfT3|) has been used. We see that S = for finite 
lattices while for N — ► 00 we obtain the established value 
6 = 3 of the standard Ising model. In the same way we 
obtain near T 



Ho(H-H s ) 
k H T> 



T 
2T> 



In 



1- 



1-8' 



N 
N-l 



M + o(s') (20) 



and we see that in both, finite and infinite lattices (with 
H a = 0) we have 6' = 1, which is different from the stan- 
dard Ising model. The finite Ehrlich-Schwoebel barrier, 
H s ^ 0, leads to 5' = 0. 

We find the spontaneous magnetization for H = near 
T c from (HSJ| 



T c /T=l + - +o(s 4 ) 



which results with (|12|) in 



(21) 



(22) 



such that we have ft = \ for the infinite-size limit and /3 = 
for the finite-size case. Near T' c we obtain analogously 
-HoH s + (±T C + \T' c )s + T' c s' = s'T + o(s' 3 ) which leads 
to 




M-- 



N-l T-Tl 



N 



X -T c 

j c 



At 1 ' 

2 c 



1 

'n 1 



N-l 



N 



k B {)T c 



•(23) 



and (3' = for finite-size lattices independent of the 
Ehrlich-Schwoebel barrier. The same exponent appears 
for infinite size since according to (|13|) the anomalous 
spins s' do not contribute to the magnetization in the 
infinite limit. 




t = T/T 



1/2 

FIG. 4: The magnetization of the Ising model (open symbols 
and dashed line) with a step as in Fig. together with the 
solution of (closed symbols and solid lines). 

Besides the divergence at T c the specific heat shows 
a second maximum at T' c = 2J'/JT C for fields fulfilling 
© as can be seen in Fig. [3J The interesting leading 
order near the critical points at H = 0, where we have 
s = y/3{T c /T - 1) from JSJ, reads 



cn/ks = 



y/3(N-l)s'T, 



3/2 



\3 



it') 

2 c) 



AN 2 (T C 



T-T c ^T~T C 



-(1-8 

9 , 1 , 
8+°^ 



)T C ') (T 



-o(VT- 



T )3/2 

T'c) (24) 



6 





a 




7 


8 


a+2/3-^y 


a+/3(H-<5) 


2D Ising (exact) 





1/8 


7/4 


15 


2 


2 


2D Ising (Weiss) 





1/2 


1 


3 


2 


2 




3/2 





1 





2.5 


1.5 


T AT 

1 c 1 V — OO 


U 


1 /o 
1/2 


1 
1 


Q 
O 


o 
z 


z 


T' c N^oo (H s =0) 


2(0) 


0(0) 


1 (1) 


0(1) 


3(1) 


2 (2) 


T' c N = oo (H s =0) 


0(0) 


0(0) 


0(0) 


0(1) 


0(0) 


0(0) 



TABLE I: Critical exponents for the two transitions in mean- 
field approximation. The value in brackets give the results 
without Ehrlich-Schwoebel barrier. 



which leads to a = 3/2 for the finite and a — for the 
infinite case. Near the other critical point T' c with s' — > 
we have the leading order 



H 



(,^H s f(l-N)[a+T' c (2T c+] r c ff o( i) 



a 2 N 2 (T-T^) 2 
+ (-r—H s jT c ) [(—H s j) 



T-T> 
{2T c +jT' c ) 2 ] 



a+2(T c -T')(2T c +jT') 2 , 1 , , , s , , 

x i en c ; +o( _ )+o(T _ r / ) (25) 



-4^fT c /k% + (T c - T'){2T C + jT' c f. It 



with a 

shows a' = 2 for the finite case and a' — for the infi- 
nite case. In the case with no Ehrlich-Schwoebel barrier 



4 



o(s' 2 ) which 
= occurs for 



(H s = 0) the specific heat becomes c' f 
shows a' = 0, and no second transition s' 
finite or infinite systems. 

One should note that the divergence of the specific 
heat appears here in the mean-field model though the 
numerical data show a mere maximum. This rounding of 
the divergence is due to the finite size of the lattice and 
well discussed, see 0- 

The results for the mean-field model are summarized 
in table U We see that for infinite lattices the presence 
of the step does not change the exponents of the Ising 
model near the normal transition T c . The finite-size ef- 
fects lead to a deviation of all exponents from the result 
without step except the exponent of the susceptibility 
which remains unchanged. Especially the specific heat 
changes from logarithmic divergence to power-law diver- 
gence. For the reported second transition the scaling in- 
equalities are fulfilled. In the infinite-size limit no second 
transition occurs. 

Please note that we compare here the mean-field crit- 
ical exponents for the new transition arising due to the 
Ehrlich-Schwoebel barrier. The exact one for the lad- 
der defects without Ehrlich-Schwoebel barrier is well 
knownS, as presented above |PJ. 



V. SUMMARY 

For Ising systems on a square lattice with a spatial 
distortion we report here that a second, low-temperature 
transition occurs besides the standard Ising phase transi- 
tion. An analytical two-spin model is capable to describe 
the main features of such a distorted finite spin system. 
The divergent heat capacity appears here due to the spa- 
tial distortion and not due to an Ising phase transition. 
Therefore, experimentally recorded signals with diver- 
gent heat capacities may not exclusively be interpreted in 
terms of phase transitions in finite systems. When sim- 
ulating the surface coverage with molecules the present 
model is able to describe the main equilibrium features^, 
thus it promises an application potential to the fabrica- 
tion of nanowires which are created near a surface step. 



APPENDIX A: CRITICAL 
EHRLICH-SCHWOEBEL BARRIER 

Here we outline the discussion of the critical Ehrlich- 
Schwoebel barrier where the second transition occurs in a 
two-spin model if condition J5| is fulfilled. A completely 
analogous discussion leads to the result for the three-spin 
model Q). 

The fictitious spin s' obeys the equation 



s' = tanh/3(# - H s 



bs') 



(Al) 



= -1 



and the second transition occurs if s = 1 and 
since only in this case the magnetization (|13f) is reduced. 
For zero temperature the tanh function shows that s' = 
±1 if and only if 



c±= H + a±b^H s 



(A2) 



Since c_ < c + we have the situation that for H s < c_ 
we have s' = 1 and for H s > c + we have s' = —1 while 
for c_ < H s < c + both solutions s' = ±1 exists. In 
this range the system will take the solution where the 
partition function becomes maximal. Since the partition 
function is proportional to cosh(c± — H s ) we have s' = ±1 
if 



H„\ % \c--H a 



(A3) 



Since we considered the range c_ < H s < c+ we obtain 
with E2l 



(A4) 

1 and s' — ±1 and where 



H + a ^ H s 

as a unique condition where s 
the second transition occurs. 
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